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Abstract 



We study the behavior of the topological defects in the inherent struc- 
tures of a two-dimensional binary Lennard-Jones system as the size dispersity 
varies. We find that topological defects arising from the particle size dis- 
persity are responsible for destabilizing the solid as follows: (i) for particle 
density p < 0.9, the solid melts through intermediate states of decreasing 
hexatic order arising from the proliferation of unbounded dislocations, (ii) 
for p > 0.9, the dislocations form grain boundaries, dividing the system into 
micro-crystallites and destroying the translational and orientational order. 

Topological defects play a crucial role in melting of a solid, especially in two dimen- 
sions. In two dimensions (2D) these defects, present in bound pairs at low temperature 
solid phase, are believed to unbind and destroy the crystalline order as temperature is raised 
causing melting M. Like temperature, size-dispersity (inhomogeneity in particle size) dis- 
favors crystalline order, and can even convert a solid to a liquid |2|,|3|]. Here we study the 
defect morphology in the inherent structures of a 2D Lennard-Jones system with a bimodal 
distribution of particle sizes. 
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We simulate N = 10 4 particles interacting with a truncated "shifted-force Lennard- 
Jones" pair-potential in 2D. We choose half the particles to be smaller than the rest, and 
define the size dispersity A to be the ratio of the difference in their sizes to the mean size. 
We start by placing the particles randomly on the sites of a triangular lattice, embedded in 
a rectangular box of edges L x and L y with aspect ratio L x /L y = \f?>/2 (to accommodate 
close-pack hexagonal structure without distortion). We apply periodic boundary condition 
and use the velocity Verlet method to integrate Newton's equation of motion. Units are set 
by choosing the mass of the particles and the Lennard- Jones (LJ) energy and length scales 
to be unity. The density p is the ratio of the area occupied by the particles to the box area. 

We equilibrate a state, defined by (p, A), at a constant temperature T = 1 using Berend- 
sen's thermostat; at this temperature the 2D solid can form for low enough dispersity. We 
run our simulation at constant energy until the temperature T, the pressure P, and the 
energy E stabilize with less than 1% fluctuation (typically for 2 x 10 5 time steps, where the 
time step is 0.01 in LJ units). We also check that the average particle displacement is at 
least a few times the average particle size. To obtain higher density states, we increase p in 
steps of 0.01 from 0.85 (liquid state) to 1.05 by gradually compressing the box, keeping the 
aspect ratio fixed, and equilibrating the system at each of these densities. 

We perform defect analysis on 100 equilibrated configurations for each of the state points 
(p, A). Clear identification of geometrical defects is difficult in simulations due to the pres- 
ence of many "virtual defects" which arise because of vibrational excitations. To overcome 
this difficulty, we analyze the inherent structure of each configuration [|J, obtained by re- 
moving the vibrational excitations — or equivalently by minimizing (locally) the potential 
energy using the conjugate gradient method. 

In order to find the defects, we construct a Voronoi cell around each particle, thereby 
uniquely defining its nearest neighbors. The ordered close-packed structure of the 2D solid 
is hexagonal, so each particle i at position r-j has rti = 6 neighbors. A defect arises when a 
particle has ^ 6, which generates a "topological charge" qi = rti — 6 ||. Defects which 
are nearest neighbors are grouped as defect clusters with total charge Q = an d total 
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dipole moment P = J2i r i<?j> where the sum is over all defects % in the cluster. Figure [l| shows 
an example of defects in an otherwise ideal triangular lattice. 

We find that the defects fall into three categories: (i) Monopoles. Clusters with Q ^ 0. 
The simplest case is a disclination — a size-one cluster (Fig. [L]a). (ii) Dipoles. Clusters with 
P ^ 0. The simplest case is a dislocation, a size- two cluster (Figs. |I|a and [TJb), composed 
of a "bound pair" of neighboring defects of opposite charge, (iii) Blobs. Clusters with 
Q — P = 0. The most common case is a quadrupole (Fig. |l]a) made of a "bound pair" of 
neighboring dislocations with oppositely-oriented dipole moments (see, e.g., 

Figure shows typical snapshots of the inherent structure at p = 0.9 for three values of 
A. The corresponding phase points (p, A) represent solid, hexatic and liquid phases 0. In 
the low-dispersity solid phase [Figs. |2|a and [|b] , defects occur mostly in the form of blobs 
(quadrupoles). Such defects, being without charge or dipole moment, cause little distortion 
in the nearby order and so are energetically inexpensive. We find that the defects in the 
inherent structure of the solid phase aggregate to form domains separated by nearly defect- 
free regions, suggesting that there is an effective attraction between defects. This effective 
attraction may arise from packing constraints since in dense packing, formation of a local 
large-amplitude defect is improbable. As A increases, more defects in the form of dipoles 
are created. 

At large defect density, free dislocations appear. Figure 0c shows their presence, in the 
hexatic phase. A dislocation destroys the long-range translational order as it introduces 
an extra half row (Fig. |l|b) that can only terminate in another dislocation with equal but 
opposite dipole moment. Translational order is destroyed over the range of separation of 
the dislocation pair Jl|. Dislocations, however, retain orientational order. We find in con- 
figurations such as Fig. that 2D translational order is lost (due to the abundance of 
dipoles). However, orientational order shows an algebraic decay — the characteristic features 
of a hexatic phase The system breaks up into crystalline patches of finite length £ t 
(the translational correlation length) which are shifted but not rotated with respect to one 
another, so the range of orientational correlation £q far exceeds On further increasing 
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A, more defects are created: many monopoles appear which destroy the orientational order 
and the system melts [Fig. |2|d] (see 

In the hexatic phase, we find a steep increase in the number of dipoles (Fig. |3|a) between 
Ash (the solid-hexatic transition value for A) and Ahl (the hexatic-liquid transition value 
for A). Also, we detect a gentler increase in the number of monopoles as the hexatic- 
liquid transition at A = Ahl is approached. KTHNY theory fl| predicts that in the liquid 
phase, the orientational correlation function C§(r) decays exponentially, C§(r) ~ e~ r ^ 6 , 
with an orientational correlation length £g that diverges as the liquid-hexatic transition is 
approached |]. Figure |3|b reveals a rapid increase of £q as the liquid-hexatic phase boundary 
is approached from the liquid side. 

Figures ^Sjo-d show typical snapshots of the inherent structure at p = 1.0, for A = 0.04, 
0.1 and 0.12. A first order solid-liquid transition is found at the value Asl ~ 0.1 0, the 
solid-liquid transition value for A. Our defect analysis shows that for A < Asl, there are free 
dislocations, and near the transition these defects line up in "strings" to form long chains of 
large-angle grain boundaries (Fig. |]c). As A — > Asl, these chains percolate, fragmenting the 
system into micro-crystallites rotated with respect to each other. Thus the grain boundaries 
simultaneously destroy both the translational and the rotational order. In the theory of 
grain-boundary-induced melting HH, the Landau free-energy expansion yields a first-order 
solid-liquid transition with the absence of any hexatic phase. We find that the formation of 
defects and the proliferation of grain boundaries occurs abruptly at Asl ~ 0.1. Within the 
resolution of our simulations, we do not see any hexatic phase 0. 

In summary, we have seen that size-dispersity induces topological defects which in turn 
destroy crystalline order, and that the mechanism of dispersity-induced melting displays sur- 
prising parallels with the mechanism proposed for the case of temperature-induced melting. 
Depending on the value of p, dispersity-induced melting can be either a first-order transition 
or a continuous transition (with an intervening hexatic phase), and the defect morphologies 
display completely different behavior in the two cases. 
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FIGURES 

FIG. 1. Identification of defects taken from our simulated system. Neutral (n« = 6) particles 
are shown as dots. Solid squares and triangles denote defects, with n% = 7 (q = +1) and n% = 5 
(q = —1) respectively. Neighboring defects form different types of defect clusters as explained 
below, (a) The large open square denotes a positive monopole. It represents a cluster of size one, 
which is the positive defect it contains, thus it is a disclination. The large open triangle denotes 
a negative monopole. It represents a cluster of size three, which consists of two negative and one 
positive defects. The arrows denote dipoles; all but the largest arrow (pointing "East") represent 
clusters of size two, and hence all but one dipole are dislocations. The largest arrow denotes a 
dipole, the vector sum of two adjacent disclinations forming a cluster of size four with a non-zero 
net dipole moment. Finally, the parallelogram contains a cluster of two n« = 5 (q = — 1) and two 
rii = 7 (q = +1) defects with zero net charge and zero net dipole moment (hence a quadrupole). (b) 
A dislocation dipole (solid arrow). Also shown is the corresponding Burger's vector (dashed arrow), 
as the vector needed to close an equilateral (in this case 3x3) contour around the dislocation. 
The thin solid line marks the extra row which ends at the dipole. The length of the dipole arrow 
represents the magnitude of the dipole P, proportional to the length of the Burger's vector. 
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FIG. 2. Defects of four different p = 0.9 configurations. The symbols used for the defects 
are same as in Figure |]. (a) A A = 0.04 configuration; this configuration is in the solid phase. 
The magnified portion includes two dislocations and a cluster of size 10 (encircled by the dashed 
line) with an overall dipole moment, (b) The inherent structure for the same A = 0.04 configu- 
ration. Most defects in (a) and (b) form quadrupoles, and some form dipoles. Notice that upon 
forming the inherent structure of a configuration, the particles with six neighbors move locally to 
form an ordered triangular lattice structure. Also many of the defect clusters, mostly low energy 
quadrupoles and some dipoles and monopoles, disappear. There remain only a few defect clus- 
ters, topological defects that cannot be annihilated by local translations performed in generating 
the inherent structure; e.g. in the magnified region, most of the defects have disappeared; there 
remains only one dislocation, the vector sum of the three previous dipoles. (c) Inherent structure 
for A = 0.06 configuration (in the hexatic phase). The arrows denote dipoles. Note that free 
dislocation dipoles appear, (d) Inherent structure for A = 0.09 configuration (in the liquid phase). 
Many disclination monopole defects, shown as open triangles and squares, appear. For clarity, 
individual defect particles and blobs (including therefore quadrupoles) are not shown in c and d. 
Dipole clusters are represented by the dipole vector for each cluster. 

FIG. 3. (a) The number of different defect clusters as a function of A along the p = 0.9 isochore. 
There is first a ten-fold increase in the number of dipoles associated with the solid-hexatic transition 
Ash- Then a large increase in the number of monopoles occurs near the hexatic-liquid transition 
Ahl- (b) Orientational correlation length as a function of A along the same p = 0.9 isochore. 
Upon approaching the hexatic phase from the liquid phase (i.e., as A — > 0.06 + ), £q diverges. 
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FIG. 4. Defects of the p = 1.0 system, (a) A A = 0.04 configuration, which is in the solid 
phase. The defects are magnified in the insets. This configuration contains only two dipoles and 
one low energy quadrupole. Note that the leftmost and rightmost defects (two dipoles) cancel, 
a feature characteristic of the solid phase, (b) The quadrupole disappears upon generating the 
inherent structure, (c) Inherent structure for a A = 0.095 configuration which is in the solid phase 
near the melting point, (d) Inherent structure for a A = 0.115 configuration, which is in the 
liquid phase. Note that the defects form lines (grain boundaries) that partition the system into 
sub-regions of solid order (micro-crystallites). Unlike the liquid in Fig. |2|d, there are few disclination 
monopoles. What destroys the solid order is the existence of the grain boundaries, while in Fig. ||d 
what destroys the solid order is the existence of both free dislocations and free disclinations. 
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